Measuring the Superfluid Fraction of an Ultracold Atomic Gas 
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We propose a method to measure the superfluid fraction of an atomic gas. The method involves the 
use of a vector potential generated by optical beams with non-zero angular momentum to simulate 
uniform rotation. The induced change in angular momentum of the atomic gas can be measured 
spectroscopically. This allows a direct determination of the superfluid fraction. 
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Central to the understanding of the physics of degen- 
erate Bose gases are the concepts of Bose-Einstein con- 
densation and superfluidity [l], Q . Bose-Einstein conden- 
sation refers to the macroscopic occupation of a single 
quantum state. Superfluidity refers to a set of fascinat- 
ing hydrodynamic phenomena, notably persistent (dissi- 
pationless) flow. 

Both phenomena admit clear quantitative definitions, 
allowing a Bose gas to be characterised by "condensate" 
and "superfluid" fractions [J, E|. These two quantities 
in general take very different values. A gas of non- 
interacting bosons at low temperature forms a Bose- 
Einstein condensate (BEC), but is not superfluid. In the 
low-temperature limit, liquid 4 He is both a BEC and a 
superfluid, but the condensate and superfluid fractions 
are markedly different, believed to be ~ 10% and 100%, 
respectively [![. In 2D systems the superfluid fraction 
can be non-zero even if the condensate fraction vanishes. 

For ultra-cold atomic Bose gases, the condensate frac- 
tion is readily measured through the mapping of occu- 
pation numbers in momentum space to real space by 
expansion imaging [4j. Characteristic signatures of su- 
perfluidity have been observed in atomic gases, notably 
dissipationless flow [H, and the formation of quan- 
tized vortices in rotating gases Q- However, there has 
been no quantitative measurement of the superfluid frac- 
tion. Such a measurement is crucial for the investigation 
of some of the most interesting properties of interact- 
ing Bose gases [8|: strong interactions can lead to con- 
densate depletion without loss of superfluid fraction; the 
Kosterlitz-Thouless phase transition in a quasi-2D ge- 
ometry is manifest by a universal jump of the superfluid 
density ||. 

In this Letter, we describe how the superfluid fraction 
of an atomic gas can be measured using a light-induced 
vector potential Our method is closely analogous to 
the classic experimental method of Andronikashvili fioj ] . 
There, liquid helium is put in contact with a rotating 
object: the normal fluid picks up non-zero angular mo- 
mentum, while the superfluid acquires no angular mo- 
mentum. A measurement of the angular momentum of 
the fluid then allows a determination of the superfluid 
fraction. Here, we consider the use of an optically in- 
duced vector potential to simulate uniform rotation Q. 



We show how spectroscopy can be used to measure the 
net change in angular momentum of the fluid, and hence 
the superfluid fraction. Recently, measurement of the 
superfluid fraction from the density profile of a rotating 
gas was also proposed [Hj]. Throughout the paper, we 
consider a gas of identical bosons, but the method can 
be extended to other situations, such as superfluidity of 
paired fermions. 

The definition of the superfluid fraction was expressed 
in a form suitable for our purposes by Leggett [13]. 
It applies to a fluid contained in a ring-shaped vessel 
with a radius R that is large compared to its transverse 
dimensions, so that the classical moment of inertia is 
I c \ = NMR 2 for N atoms of mass M. We start by adopt- 
ing this assumption of geometry, but this will be relaxed 
at the end of the paper. The walls of the vessel are taken 
to rotate with angular velocity u>, and the fluid allowed 
to come to thermal equilibrium. Under these conditions, 
the superfluid fraction is [l2[ 
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where (L) is the average angular momentum of the fluid. 
A normal fluid will come to rest in the frame rotating 
with the walls, so that (L) = I c \lo and p s /p = 0. A 
(perfect) superfluid is unaffected by the rotating walls, 
so (L) = and p s /p = 1. 

When the fluid is in equilibrium with the rotating walls 
it is described by the equilibrium density matrix defined 
by the Hamiltonian in the rotating frame [131 ] 

H Iot =H-L-u. (2) 

Here, H and L are the Hamiltonian and the total angular 
momentum in the laboratory frame. We shall show how a 
Hamiltonian of the form ([2]) can be achieved for an atomic 
gas, and how the resulting average angular momentum 
(L) can be measured so that (fTJ) can be directly applied. 

In the ring geometry, the kinetic energy in ([2]) can be 
written 
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where p± is the momentum in directions perpendicular 
to the azimuthal direction, and I is the angular momen- 
tum in units of H (therefore quantized to integer values). 
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The rotation shifts the energy minimum in the angular 
momentum to 
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This shift can be viewed as an azimuthal vector potential 
corresponding to a non-zero flux threading the ring. 

A shift in the dispersion relation can be achieved by the 
use of two-photon Raman transitions to imprint vector 
potentials [8(. This was recently implemented 14 1 us- 
ing two counter-propagating laser beams to couple states 
m = — 1, 0, 1 of the F = 1 hyperfme levels of 87 Rb. The 
two-photon processes lead to a linear vector potential di- 
rected along the axis of the lasers [3, EH ■ 

To generate an azimuthal vector potential, we consider 
two Laguerre-Gauss (L-G) beams [16| with different or- 
bital angular momenta, co-propagating in the direction 
perpendicular to the toroidal trap. In this way, a two- 
photon transition imparts negligible linear momentum 
to the atoms, but a non-zero angular momentum, ±A£, 
where A£ is the difference in the orbital angular momenta 
of the two beams. For a 3-level system 14, lj| this leads 
to an effective Hamiltonian 
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which is a matrix in the hyperfme states m = —1,0, +1. 
flu is the two-photon Rabi frequency, S = is 
the detuning of the lasers from the Raman resonance set 
by the Zeeman effect of a (uniform) magnetic field B, 
and e accounts for the quadratic Zeeman effect. 

The energy eigenvalues for the angular motion are il- 
lustrated in Fig. [TJ The lowest energy band has a mini- 
mum at a non-zero angular momentum £* . To derive an- 
alytic expressions for the energies and wavefunctions in 
this band, we consider the limit of large Rabi frequency 
and develop a perturbation theory in 1/£Ir. We param- 
eterise the resulting single-particle energy as 
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The energy En is a global shift that depends on all pa- 
rameters of the optical field. A slight increase in the 
effective mass for the azimuthal motion is given by 
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The most important effect for our purposes is the shift 
in the minimum of the dispersion curve to 
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FIG. 1: Energy levels for the angular motion of an 
atom under the influence of two-photon Raman coupling 
via beams with orbital angular momentum difference At. 
The atoms move in a trap of radius R, the detuning is 
S = 0.5h(A£) 2 /MR 2 , e = 0, and Rabi frequencies of Q R = 
(dotted line) and Q R = 2h(A£) 2 /MR 2 (solid lines) are shown. 
The lowest band has its minimum displaced to a nonzero an- 
gular momentum, £*, equivalent to the effect of an azimuthal 
vector potential. (Note that the smooth curves should be 
viewed as interpolations between the allowed integer values 
of £.) 



This can be viewed as the introduction of an effective 
vector potential. In order to clarify the physical conse- 
quences of I*, note that, for an atom with angular mo- 
mentum £, the angular velocity when the light is on is 
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Thus, for given £, the most significant effect of the light 
is to cause a constant shift in the angular velocity by 
h£*/M*R 2 . In analogy with f3j|4]), it is as if the optical 
field causes the laboratory frame to behave as a frame of 
reference that is rotating with angular frequency 
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The lowest band (JSJ) plays the role of H Iot ((2), with 
an effective rotation rate © that can be tuned by the 
parameters of the optical field. Provided the splitting of 
the bands is large compared to the chemical potential, 
^ M/% an " atoms are restricted to this lowest band. 
Then, if the position of its minimum |[7J) is varied suffi- 
ciently slowly that the fluid has time to come to equilib- 
rium for the new £* , a clear distinction appears between 
normal and supcrfluid components: the normal fluid will 
relax, and pick up a non-zero angular momentum; the 
superfluid will not relax, but will retain vanishing aver- 
age angular momentum. This distinction provides the 
definition of the superfluid fraction |T]). 

This behaviour of the atomic gas in a toroidal trap 
with an azimuthal vector potential is in marked contrast 
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to the case of a linear vector potential [14[. There, both 
superfLuid and normal components must come to rest in 
the laboratory frame, consistent with the observations in 
Ref. [3]. Here, when the normal fluid comes to equi- 
librium at a non-zero £*, it is at rest in the laboratory 
frame. However, the superfluid is rotating, as follows 
from ([8]) . The azimuthal vector potential causes a steady 
superfluid flow around the ring-shaped trap. Note that 
for the normal fluid to come to equilibrium with the new 
£*, it must change its angular momentum. Therefore, the 
trap must not be perfectly rotationally symmetric. (No 
fluid can come to equilibrium in a rotating container if 
the walls of the container are perfectly smooth.) It is an 
important practical feature that there is no requirement 
for the trap to have perfect rotational symmetry. 

In the Andronikashvili experiment [1 01 ] a torsional os- 
cillator is used to measure the moment of inertia of the 
fluid coupled to the oscillator, /. A frequency shift arises 
from the fluid's contribution to the energy, which in the 
rotating frame is (H rot ) = — \Ilo 2 . Here, when the light is 
on there is a contribution to the energy of — ^Iu> 2 s . One 
can envisage various ways in which an oscillator can mod- 
ulate w c f ([9]) and therefore experience a frequency shift 
related to /. For example, given that oj c s « f a B, the 
coil that generates the Zeeman field B experiences the 
moment of inertia as a reduction in its inductance; this 
could appear as a shift in the resonant frequency of an 
electrical circuit containing the coil. However, the total 
energy \l(TiA£/M* R 2 ) 2 is very small (~ 0.1/xeV) mak- 
ing the signal small compared to typical sensitivities of 
current micro-mechanical or electrical oscillators. 

A key element of our proposal is that, with the 
above light-induced vector potential, one can use spec- 
troscopic methods to determine the average angular mo- 
mentum (L). The wavefunction in the lowest band 
is a linear superposition of the three hyperfine levels 
= J2 m =-i,o,i ^m| m ) wrtn amplitudes {ip m } which 
vary with I. A perturbative analysis shows that there 
are equal and opposite corrections to (i/'iil 2 which de- 
pend linearly on £. Thus, (L) can be obtained from a 
measurement of the difference in the number of particles 
in the states m = ±1. Using the Taylor expansion 

| 2 - |V^i| 2 = Ap + Ap'£ + O(£ 2 ) (10) 

we can write 

(L) _ Z e (n e )£ _ Ap-Ap 
where we define the fractional imbalance 

A jV-i-JVi = E/("<) [hM 2 -hAi| 2 ] (12) 

with (ni) the average number of particles with angular 
momentum £. Inserting (|11[) in JT]) one finds 

(3^) +<*•/«»«>■ 



where the corrections arise from the approximation (|10p , 
which is accurate provided the atoms are in the parabolic 
region of the lowest band. 

In the limit of large Q R , the wavefunction of the lowest 
band is (-0-1,^0,^0 = 1/2(1, -\/2, 1) for all £. Com- 
puting perturbative corrections to order 1/fifj in l^ml 2 , 
we find A Po ~ (S/n%) [V2Q R - h(A£) 2 /(2MR 2 ) - 2e] + 

o(i/n%), and a p ' = V2hA£/(MR 2 n B ) + o{i/n 2 R ). 

Combining this with ([7]) we obtain 

The limit u> — » in ([T]), replaced here by to e s oc S — > 0, is 
discussed further below. 

Eqns. (|131 1 14|) show how a spectroscopic measurement 
of the populations N m can lead to a direct measurement 
of the superfluid fraction. That there is a connection 
between these quantities is a central result of this paper. 

To distinguish a normal fluid from superfluid, the frac- 
tional population difference Ap (|12[) must be measured 
with an absolute accuracy of order 

2MAI) 2 S 

mr 2 n 2 R ■ [Lb) 

This expression was derived for 6/CIr <C 1. In Fig. [2] we 
show the expected fractional change in occupation for 
a normal fluid, with angular momentum centred on £*, 
computed for arbitrary S/Q R . This is shown for param- 
eters which for 23 Na would correspond to R = 10 pm, 
fl R ~ 2n x 4.4 kHz, and A£ = 10 (two beams of or- 
bital angular momentum 5). In this case, Ap must be 
measured to an absolute accuracy of about 3%. The 
required relative accuracy to distinguish a normal fluid 
(Ap — Apo 7^ 0) from a superfluid (Ap = Ap ) is 
(Ap — Apo)/Apo, and is about 10% in the linear regime 
8/VLr < 0.25 in Fig.O This relatively small signal poses 
a moderate experimental challenge. It is important to 
stress that it relies only on the measurement of fractional 
occupations of different states. It is therefore insensitive 
to systematic uncertainties in the absolute atom number 
determination, and statistical errors can be reduced by 
averaging over many shots. 

Above we used A£ = 10 as a currently realistic experi- 
mental value, but the signal can be increased by increas- 
ing A£. Since we require 11 r > [i/ft, for atoms to be in the 
parabolic region of the lowest band, it is instructive to 
write ([IS]) in the form (n/hn R )(A6/n R )(A£^/ R) 2 , where 
£ is the healing length. All three terms in this expression 
can in principle be close to unity. One limitation to A£ 
arises from the fact that the expression (|14p applies in 
the limit Lj c g cx 6 — > 0, so that the imposed rotation is 
sufficiently small that the superfluid does not change its 
angular momentum. This will fail if the resulting veloc- 
ity of the superfluid, Rui e s = h£* / (M*R), becomes larger 
than the superfluid critical velocity, which is ~ K/(M£). 
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FIG. 2: Angular momentum £* at the bottom of the band 
(dashed line) and change in particle imbalance Ap—Apo (solid 
line) as a function of &/£Ir for a normal fluid (i.e. centred on 
t) for Q, R = 1000ft/ MR 2 , At = 10, e = 0. This illustrates 
the precision required to distinguish a normal fluid (here Ap — 
Apo ~ 3%) from a perfect superfluid (Ap — Apo = 0). 

If this condition is violated, the superfluid will relax (vor- 
tices will enter the system) and acquire non-zero angu- 
lar momentum. The condition for stability of the super- 
fluid flow requires A£ to be less than ~ (SIr/5)R/£ for 
S <C 0, and ~ R/£ otherwise. (Note that typically 
£ < 0.5 /urn.) As (fl~5|) suggests, if A£ is limited by practi- 
cal reasons, the signal will generally be larger for lighter 
species, which for typical experimental parameters have 
lower density, and hence larger £. 

Experimentally, a large spectroscopic signal which 
qualitatively distinguishes normal and superfluid compo- 
nents can be observed by a "projective measurement": 
Suppose that the system is in equilibrium in the parabolic 
band with Qr ^> » 5. If we then reduce Qr to ~ S 
on a time scale short compared to the relaxation time but 
long compared to 1/VLr, the superfluid and the normal 
fluids will remain centred at I = and I = £*, respec- 
tively, but the difference in their spin composition will be 
greatly enhanced (see Fig. [I}. However the quantitative 
extraction of the superfluid fraction would in this case 
require further analysis. 

The ring-shaped trap discussed above is the case clos- 
est to the theoretical discussions of the superfluid frac- 
tion and the simplest to present. However, our 
method also applies to a quasi-2D or 3D gas, provided 
the optical fields are such that the atoms always remain 
in the lowest energy band. As an illustration we con- 
sider a scenario in which two hyperfine levels [15j, la- 
belled | and J., are coupled by the L-G beams propagating 
along z. In the (ipi,i/>\.) basis, we parameterise the low- 
est energy dressed eigenstate as (e lx sin(0/2), cos(0/2)), 
where 9(r) and x( r ) depend on the local optical field. 
For beams with angular momentum difference At, we 
take x{ r ) = A£(j), where <j) is the azimuthal angle around 
the z axis. Assuming optical fields such that 9(r) = 



ar <C 1, where r is the radial distance in the xy plane, 
the lowest energy state experiences an effective vector 
potential that simulates uniform rotation ||. The to- 
tal number of flux quan ta inside r is A£(ar /2) 2 and 
^eff = (h/4M) A£a 2 [V^. Computing the leading per- 
turbative corrections to the state, as in (flT))) . one finds 
IV'tI 2 ~ IV'il 2 - &Po + Ap'r(p^(r)>, where (p^(r)) is the 
local azimuthal momentum density. Thus, the correc- 
tion to iV'tl 2 — IV'il 2 depends linearly on the local angular 
momentum density. The spectroscopic measurement of 
(JV-f — Ni)/N, integrated over the sample, therefore pro- 
vides a measure of the angular momentum per particle 
(L) /N. Comparison of (L) with I c \uj e s allows the deter- 
mination of the normal and superfluid fractions. 

In summary, we have proposed a method to measure 
the superfluid fraction of an ultracold atomic gas. It 
combines the use of optical beams with non-zero angu- 
lar momentum to simulate rotation, with a spectroscopic 
readout of angular momentum. Our observation that 
light-induced vector potentials create a direct connection 
between the formal definition of superfluid density and 
the spin composition of a gas is very general, and we ex- 
pect it to be applicable to other experimental scenarios. 
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